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1.  INTRODUCTION 


In  the  process  of  evaluating  various  turbulent  mixing  models,  the  question  kept 
arising  of  how  the  turbulent  flow  differed  from  the  laminar  flow.  A  comparison  of  two 
different  turbulence  kinetic  energy  models  is  given  in  Walker  [4]  and,  in  addition,  the  use  of  a 
laminar  viscosity  model  for  comparison  with  the  turbulence  models  is  discussed.  This  was 
accomplished  for  a  reacting  shear  layer  for  which  experimental  results  had  been  obtained. 
Therefore,  the  techniques  and  results  presented  herein  were  utilized  to  make  qualitative 
comparisons. 

Additionally  the  techniques  presented  in  this  work  are  directly  useful  in  the  chemical 
laser  program  because  of  the  fine  scale  mixing  in  the  nozzles  and  the  low  pressure  operation  in 
the  laser  cavity.  The  coding  presented  herein  is  applicable  only  for  a  N2,  O2,  NO,  NO2,  and  O3 
system  but  minor  coding  changes  will  make  these  results  applicable  for  any  gas  at  low 
pressure. 

II.  CUBIC  SPLINE  INTERPOLATION 

In  order  to  determine  the  laminar  viscosity  of  the  various  gases  involved  as  a  function 
of  temperature,  it  is  necessary  to  interpolation  Table  I  for  temperatures  not  given.  This  could 
be  done  utilizing  a  simple  linear  interpolation  scheme.  However,  a  more  accurate  scheme  which 
does  not  take  an  excessive  amount  of  computational  time  is  cubic  spline  interpolation.  This  is 
a  piecewise  cubic  interpolation  scheme  that  matches  the  function  and  its  slope  at  each  of  the 
known  points  given  in  Table  1.  Cubic  spline  interpolation  is  currently  very  popular  and  will  be 
utilized  here. 

The  only  problem  with  this  scheme  occurs  at  the  end  points  of  the  interval.  The 
requirement  that  the  slopes  be  matched  at  this  point  between  the  piecewise  cubic  sections 
cannot  be  met  since  no  other  point  is  available  from  which  to  calculate  the  slope.  Therefore  a 
special  treatment  for  the  end  points  was  necessary  and  the  use  of  a  Lagrangian  polynomial  was 
chosen, 

The  Lagrangian  polynomial  is  given  by 

n 

p(x)  *=  Z  f(x.)  1.  (x)  (•) 

k»o 
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PREDICTION  OF  TRANSPORT  PROPERTIES  OF  GASES  AT  LOW  DENSITIES  [s] 
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Now  when  these  equations  are  used  to  calculate  the  left  end  point  slope 


Becomes  X2 
Xj^  Becomes  x^ 

and  when  they  are  used  to  calculate  the  right  end  point  slope 


x^  Becomes  x,  , 
o  k-3 

Xj^  Becomes  Xj^_2 
as  shown  in  Figure  I. 


The  special  end  point  treatment  was  incorporated  with  a  cubic  spline  interpolation 
routine  developed  by  Conte  and  de  Boor  [  1  ].  The  program  listing  is  given  in  Appendix  A.  As  a 
test  for  this  generalized  routine,  the  viscosity  function  was  interpolated  at  several  points 
including  input  points  to  determine  how  well  the  cubic  spline  interpolation  program  was 
working.  Table  I  gives  the  known  input  data  for  as  a  function  of  xT/e. 


The  results  of  the  use  of  the  cubic  spline  interpolation  program  are  given  in  Table  2. 
Note  that  the  interpolation  scheme  gives  very  smooth  results.  Also  note  that  the  input  data  is 
under  no  restriction  as  to  the  regularity  of  the  input  function  interval  i.e.,  A(»cT  /e)  varies  from 
0.05-5.0. 

III.  MOLECULAR  VISCOSITY  OF  NO2  AND  O3 

The  reacting  shear  layer  that  calculations  were  made  for  contained  the  following 
reaction: 

NO  +  O3  -  NO:  +  O2 

Since  the  mixed  stream  includes  a  combination  of  these  gases,  it  is  necessary  to  know  the 
laminar  viscosities  of  each  of  the  component  gases  as  a  function  of  temperature.  These  have 
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la  .1AToOAAOF*0 
la  .lAPnOOOQF.O 
la  .)aoOAOOOf*A. 
la  .140nAAOOf*af 
la  .IMIPOPAOC.AI 
la  .1A?OAO00E*AI 
la  .1O3OA00aE4Al 
la  .  .IO4QAAA0F401 
la  ,10MO<IAOOE«A1 
la  .|06AORaoe*Ol 
|a  ,10TOAOOAI«0 
la  .10ROOAOOE*A 
la  .1O«0000AI*P) 
la  .170AAOOOF*A1 
la  ,I7IOAOOOE*01 
la  .ITPOOAOOFnOl 
la  .173AP0A3P.0J 
la  . 17400, )OQE*a 
l«  . j7aoOOPOF»Ol 
la  .ITnOQOOOff^O 
la  .ITTOOOOOFaO 
la  .ITMOOOOOFnA 
la  .iTVAROnOF*!)) 
la  .1POPOO0OF*0t 
la  .IPIOOOfOEnO 
la  ,1OP0OO0nE*.0 
la  .10300000F*0 
la  ,1P400000E*0 
>a  .1R<tOO(IOOF*01 
la  .|POOOO00F*0 
la  .10700000^*01 
|a  ,1AOoaoOOE*C , 
la  ,ia9000onF*P 
la  .IM000O00C*Ol 
la  .iPioononF*oi 


F  Ilia 
Fllia 
FUta 
F  |l»a 
Fllia 
F  111* 
Fllia 
F  Ilia 
Fllia 
F|l»a 
FlUa 
F  |l)a 
F|l|a 
F  |l)a 
F  |*)a 
F|l)a 
Fllia 
F|l)a 
F|l)a 
F  |l|a 
Fllia 
F|l)a 
F|l)a 
Fllia 
F|l|a 
FHia 
Fllia 
Fllia 
F|l|a 
F|l)a 
F  liila 
Fllia 
F|»»a 
F|l|a 
Fllia 
F|l|a 
F|l)a 
Fllia 
F  |»)a 
F|l)a 
Fllia 

7!!!: 

f!;;: 

FHia 

Fllia 

F|l)« 

F|l|» 

F|»)« 

Fllia 

F(l)a 

Fill- 

Fllia 

Fllia 

Fllia 
F  l»ia 
Fill* 
F  |l|a 
Fllia 
F|l|a 
F  ll)a 
F|l»a 
F  ll)a 
Fllia 
F|l»a 
Fllia 
Fllia 
F|l)a 
F|l|a 
Fllia 
F|l)a 

Fin* 

Fill* 

F|l»a 
F  Ilia 
Fllia 
Fllia 
F|l)» 
Fllia 


.|S0T3T<>3F*e: 

.i«oohp??f*o: 

.14M44Msme*o; 

.I4MM?0S8F*e' 

.i40PooAo^*o: 

.147500071*0 

.i4(>9op?or*o! 

. 1 463004 a|*0 
.1457*7061*0 
.|4SP0OA0E*e 
.I4461940r*e 
.14404090^*0 
.|43404H7F*0 

.lA^MSaiTE*© 

.14?400A0|*0. 

.14l00054r*0i 

.141374?7C*0* 

.140HT773e*0 

.14030746C*0 

:}?r4?5J^:8 

.1309?7)?F*0 
.l3«4450nF*0 
.137960591. A 
.137S0P0AE*0 
.13704096E*0 
.136S916OF*0 
.l361SiP6E*0 
.13S7213HF*0 
.l3630OAAr*0 
.134807341*0. 
.13440PP9F*0| 
134A03S7E*e 
|336H990E*0 
‘333000AE*0 

iEioI 

31771 63E*0 
3)400A0F*t 
3taS3O9C*0 
i067P4W40 

PPOipPM*0 
P9P40T6E»e 

?i»S97?9E*0 
?e55443E*0 
PHPE|06E*e 
E79API>a|*A 
P759POOF.0 
.j?7P9470F*0 
.|?7A0nS4C*A 
.I?O704?AF*0 
.|P6400()AE*0, 
.1P„0041  tE*0:l 
.1?970039E*0 
.1?543441F*# 
.i?5m55e*o 
.iP4i>oonoF.*o 
.1?45eiPHF*0 
.1?4?I407F*0 


.ll2U«7?F*0i 

.l?260046F*e| 

.1??3M*6E*0 

.1P2100POF*0 

,121*S477E*0 

iiWimi 

.1?090nOA|*Ol 

.1?A650JPE*0* 

.)?041400r*i 

:}?Sii5S8f:8 

:ir3§§r.8 
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TABLE  2.  (CONTINUED) 


PT.MO.a 

163 

>fe 

PT.Nft.* 

164 

Vs 

PT.no.* 

16b 

JTs 

PT.nn.a 

166 

Va 

PT.no.* 

167 

ys 

PT.nf).* 

16B 

Vs 

PT.Wf).* 

169 

Xs 

OT.wn.a 

1  70 

Xa 

PT.WO.* 

171 

Us 

PT.no.* 

17? 

Xa 

PT.WO. a 

173 

Xa 

PT.WO, a 

174 

Xa 

PT.no.* 

179 

Xa 

PT.WO. a 

176 

Xa 

PT.WO.* 

177 

Xa 

PT.no.* 

17fl 

Xa 

PT.NO.* 

179 

Xa 

PT.WO.* 

1«0 

Xa 

PT.WO.* 

1B1 

-  Xa 

PT.no.* 

)R? 

Va 

PT.no. 3 

1K3 

Xa 

PT.no, a 

1«4 

Xa 

PT.WO. a 

l«b 

Xa 

PT.WO. a 

1B6 

X« 

PT.no. a 

1P7 

x« 

PT.WO. a 

IBB 

Xa 

PT.no. a 

1B9 

Xa 

PT.WO. a 

190 

Xa 

PT  .WO.a 

191 

Xa 

PT.wO.a 

19? 

Xa 

PT.WO. a 

193 

Xa 

PT.WO. a 

194 

Xa 

PT.WO. a 

195 

Xa 

PT.WO. a 

196 

Xa 

PT.NO. a 

197 

Xa 

.PT.no. a 

19fi 

Xa 

PT.nn.a 

199 

Xa 

PT.WO.a 

?0Q 

7* 

PT.nn.a 

?ni 

Xa 

PT.wO.a 

?0? 

X* 

PT.wO.a 

?03 

Xa 

PT.no, a 

?04 

Xa 

PT.wO.a 

?05 

Xa 

PTiWO^a 

?06 

Xa 

PT.no. a 

?07 

Xa 

PT.WO  .'8. 

?Ofl 

Xa 

PT.no. a 

?09 

X* 

PT.WO.* 

PIO 

Xa 

PT.WO.*  , 

Pll 

X* 

PT.WO.* 

PI? 

X* 

PT.wO.a 

pn 

Xa 

PT.NO.» 

PT.N01.S 

PT.wn.a 

PT.HO.a 

PT.«n»* 

PT X 
PT,Kn,s 
PT.wn.x 
PT.no.* 
PT.MO.s 
PT.MO.i 
PT.nn.a 
PT.PO.x 

PT.no.* 
PT.VO.* 
PT.NO,* 
PT.NO.a 
PT.no.* 
PT.no,* 
PT.NO.* 
PT.WO,* 
PT.NO.* 
PT.WO.* 
PT.NO.* 
PT.VO.- 
PT .no, a 
PT.WO.* 
PT.NO.* 
'PT.WO.* 
PT.WO.* 
PT.no,* 
PT.NO,* 


.ip/ppoooc.ni 
.i«*?onoonp.oi 
.  l<>*noooof:.0) 
.inSf.ftnooF.rti 


.igvnnonoF.O) 

.PoononOftF.ni 

.PftlOftOftftF.OJ 

.PopooftOftp.ftl 


.pnnoonooF.fti 

.pftooooftoF.ni 

.P07ftOOOOF.ni 

.pftoooonoF.oi 


*  .?! POOftOOF.Ol 
=  .PI 30ftonnF.fti 
a  .P14ftOflOOF.01 

*  .PI‘iftOOOOE.01 

*  .?l#>ft0000E.01 
=  .pjTonnooe.oi 

*  .PIHOOnOOF.ftl 

*  .PISOOOftOE.ftl 

*  .ppooooooe.oi 

=  .??! ftftftOOF.Ol 
=  .PpPOOftOOE.Ol 
=  .P?3OOft0OF.Ol 
a  ,.PP4OOn00E.OI 

*  .PPSOOOOOF.Ol 

=  .??#>noftnoF..oi 

*  .PPTftftftOOE.Ol 
=  .PPOOOOOOE.OI 

*  .P?‘»OOOOOF..O\ 

*  .P.30flnoooE*n) 
»  .P3100000E.01 
a  ,?3?OOnfl0F*Ol 
a  *.?3.300000E.O) 
a  .PPOOOOOOE.OI 

*  .?3POOOOOE*ftl 

*  .?3Oft0000E.01 

*  ,?3700000E.Ol 


X*  .P4100000E.01 

)f*  .p«?nooooE*oi 

X*  .?*3nooooE.oi 

X*  .P44ooftOOE.ni 
Xa  .P4POOnO0P.O1 
Xa  .P4O00000E*ni 
X*  .?47noftooe.oi 
X*  .P4PflnftOoF.oi 
Xa  .P49pon00E.Ol 
x*  .PsoflnoooE.oi 
X*  .PSIftooOOF.OI 
X*  .P9?oooooE.ni 
Xa  .?930ftOOOE.ni 

. Xa  .p940ononF.ni 

*=  ,P990ftno6F.O) 
X*  .PsnoftnooF.O) 
Xa  .p970ftno-nE»ni 
*=  .papnonooF.fti 
X*  .pavnooooF.oi 
X*  .pnoonoooE.oi 
Xa  .poinonnoF.oi 

Xa  .P4pnnoooF»oi 

Xa  .P03onoonE*fti 

Xa  .PA4nOnOOE.01 
Xa  .?f>90000nE*Q\ 
X*  .P4f>oOftOOE.Ol 
Xa  .PoxnarnnoE.fli 
Xa  .POMoOnOOE.OJ 
X*  .phvoonnoF.ni 
Xa  .PTonnnooE.fti 
X*  .PTionnooE.oi 
X*  .prpoonooF.oi 
X*  .P73oooooE*ni 
X*  .P74nnoooF»oi 
X*  ,?7Sft0000E.01 


E<X>a  .UOPSPOOE.ftJ 
E(x>=  .1 1S03494F.0  1 
FIX)*  .1 1HH1K17E.01 
FIX)*  .IIHhOftOOF.Or 
FIX)*  .  j  )«37H(«7F.ft 
F l»)a  .IlHlSnlHF.O  1 
f |X)a  .1 17934nSF.01 
F|X)=  .l)7714hlE.01 
F(X)a  .1 ITSOOftOE.Ol 
FIX)*  ,n7?9105E.0 
FIX)*  .1 170»9«7£.01 
FIX)*  .11OB93P4E.0 
F|X)a  .llOTOlleF.O 
F|X)a  .114S1PPPE.0 
F|X)a  .)  lr>3P7SSF*n 
FIX)*  .1  lM444nE.O 
FIX)*  .115VO?4?F.O 
FIX)*  .)  l!>7h)4?F.O 
FIX)*  .iisoonnoE.o 

FIX)*  .nS41777E.O 

FIX)*  . ) J&P3493E.0 

FIX)*  .llSOSlOOE.Ol 
FIX)*  .1 14Rft9ft9F.01 
FIX)*  .1  ]4f>049PF»or 
FIX)*  .)]4S0S79E*0I 
FIX)*  .) 143P013E.01 
FIX)*  .n414P33E.01 
FIX)*  .n397?OpF.ftl 
FIX)*  .1 i38onnnF.oi 
FIX)*  .1  13f.3niBE.0l 
F  |X)s  .1 134f,3PPF.01 
F|X)*  .1 13?9Bf>f.F.01 
F|X)»  .1 1313BBBF.01 
FIX)*  .nP97494F*01 
FIX)*  .1 1PB1BH7F.01 
F|X)a  .1  lP6f.P40F.01 
F  IX)  a  ,1  12507?f.E.01 
F|X)a  .1  lP3S3PlP.nlt 
FIX)*  .1 iPPonnpE.oi 
FlX)a  .nPD4741F.01 
FIX)*  .1 llB9S4nF.01 
FIX)*  .1  U74397E.pi 
F (X)a  .11159311E.01 
FIX)*  .1  n44PB3F.01 
F(X)«  .lilP9313F.oi 


F(X)«  .inP9313F.0l 

FIX)*  .niiosvoEroi 
F|X)=  .1  lft99S43F.l11 
F|X)«  .1 10g4743F.01 
FIX)*  .nOTOOOOE.Ol 
FIX)*  .1 105531BE.01 
F|X)a  .1 lOOOTlBF.PI 
F  (X)a  .11026PP7F..01 
F  IX)*  .MOIIRXOE.OI 
FIX)*  .)e997F73F.01 
FIX)*  .109H3663F.01 
FIX)*  .)DSlf>9B64F.O} 
FIX)*  .1fl956304F.01 
FIX)*  .109430n7E.OI 
FIX)*  .10930000F.01 
FIX)*  .)0917?97E.01 
F|X)a  .)0904f<67E.01 
F|X)a  .)0H9?46SF.0) 
FIX)*  .10BB0F4MF.01 
FIX)*  .)nB6B773F.01 
FIX)*  .)n8569«.6F.o) 
FIX)*  .1 nM4SP74F.01 
FIX)*  .10H33S63F.01 


FIX)*  .10«21H?0E.ni 
FIX)*  .lOBlOnftOF.Ol 
FIX)*  .1 nT9B073F.oi 
F  IX)*'  .)07B60SSF»o1 
FIX)*  .1  0773973E’.01 
FIX)*  .)O761BS7F.0i 


FIX)*  .)0761BS7F.O , 
FIX)*  .in749T33F.01 
FIX)*  .)0737631E.0l 
Fl*)»  .10725ST9E.0 

Fix)»  .in69onnoE.o 
FIX)*  .1  0f.7n4p]E.0 
F|X)*  .1066699BF«a 
FIX)*  .10F.557I3E.ni 
FIX)*  .10644S69F.n 
FIX)*  .in633543E.ni 
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TABLE  2.  (CONTINUED) 


PT ,8n,« 

^47 

Xe 

PT.K,n,3 

2*8 

ys 

“T.mO.= 

740 

ys 

PT.NO.s 

750 

Xc 

PT.NO.x 

2S1 

ya 

PT.»in,= 

252 

Xs 

PT.Mfl.* 

253 

X* 

PT.Mrt.r 

25* 

ya 

PT.80.S 

255 

ys 

PT.Nn.» 

258 

Hm 

PT .MO.s 

257 

ys 

PT.MO.s 

258 

Xs 

PT .ND.s 

259 

X* 

PT.NO.o 

280 

Xs 

PT.MO.b 

281 

Xa 

PT.NO.b 

282 

Xs 

PT.NO.« 

283 

X« 

PT..NP.4 

28* 

y* 

PT.PO.o 

285 

ya 

PT.Un.a 

286 

ya 

PT.Nrt.* 

287 

Xa 

PT.NO.x 

288 

ys 

PT.wn.« 

289 

ya 

PT.MO.* 

270 

xa 

PT.NO.b 

271 

vs 

PT.on,* 

272 

ya 

PT,NO.« 

273 

ya 

pt.no.b 

27* 

ya 

PT.NO.s 

275 

y« 

PT.*JO.= 

278 

Va 

PT.Urt.* 

277 

ya 

PT.mo.x 

278 

ya 

PT .HP,* 

279 

Va 

PT.PP.s 

280 

vs 

pt.mo.* 

281 

ya 

PT.MO,* 

282 

xa 

PT,NO,« 

283 

xa 

PT.Nn.x 

28* 

xa 

PT.MO,* 

285 

xa 

PT.NP,« 

288 

xa 

PT.H(1.» 

287 

xa 

“T.nn,* 

288 

xa 

PT.MO.* 

289 

xa 

PT.MO.* 

200 

xa 

PT,Wfl,« 

291 

xa 

PT.Nrt.* 

292 

xa 

PT.Nft.* 

293 

xa 

PT.HP.m 

29* 

ya 

OT.MP.o 

295 

ya 

PT.MO.o 

298 

xa 

PT.NP,» 

297 

Xa 

•PT,NO,« 

298 

Xa 

PT.MO.* 

299 

Xa 

PT.MO.* 

300 

Xa 

PT.MO.« 

301 

Xa 

PT.Nfl,* 

302 

Xa 

PT.NO.* 

303- 

Xa 

pt.np.* 

30* 

xa 

PT.90.» 

305 

xa 

pt.no.« 

308 

xa 

PT.wn,» 

307 

Xa 

PT.WO.a 

308 

xa 

PT,*H1,« 

309 

Xa 

PT.NO.o 

310 

xa 

PT.NO.o 

311 

Xa 

PT .HO.* 

312 

Xa 

PT.NO.a 

313 

Xa 

OT.WP,» 

31* 

Xa 

PT.MO.s 

315 

Xa 

PT.Hfl.* 

318 

Xa 

PT.«IO,« 

317 

ya 

PT.Mn.B 

318 

x» 

PT.r»o.» 

319 

Xa 

PT,*(0,« 

320 

Xa 

PT.Fin,* 

3?1 

Xa 

pt.ho** 

32? 

xa 

PT.KO.o 

323 

Xa 

PT.wn.o 

32* 

Xa 

PT.HO.m 

325 

.Va 

PT<wn,« 

328 

Xa 

PT.MO.a 

327 

Xa 

pT.N<».« 

329 

ya 

PT.MO.o 

.329 

Xa 

PT.Nfl.o 

23® 

Xa 

PT.W0.« 

331 

Xa 

,?77PnnonF ♦« 
,?7HnonooF*r> 
,77vnonoPK*P 
,7Pnpoft00F*n 
.?»<10flnOOE*0 
,7f>?oonPOE*ft 
.7P.70P000E*P 
.?p*nrftnoF*o 
.7P?>fiPonoF.O 
,?PFnOPOOe*P 
.7fl7OOP0OF*O 
,?phonf>ooF*o 
.7«9rOrtOOF*0 
,?<»nonnooF*o 
,7«))Pftp0OF*0 
.?o?P0ft00E*0 
.?‘»30nft00E*0 
.7<)*onnooF*o 
.pOSOOOOOE^O 
,?9PnnnoPF*o 
.?‘»7n0POOF'*0 
,?«»pnoAOoe*o 
,?Q«)00P00e»0 
.3onpnAnoC«0 
,3P1 pnnooF*n 

.io?onnooE»0 

,3fl3firtP00F»0 
.304PnoooF»o 
,3nsnopooF*n 
.30FpnpooF*o 
.30700000F»n 
.IOHOOOOOF*© 
.  3(»<iopno«F*o 
.31 0(>nftOOF*o 
.3UOOOOOE*0 
.31?OOPOOF*0 
.31 3PftnooF«o 
,3\40PnQ<lF*0 
.3l^nnnoo€*o 
.3iFnoPOOE*n 
,317flftP00F*0 
.31HPn000f*0 
.SIPOOPOOF*© 
.37n«onooF*o 
.37lO0n00F*0 
.3??POOOOF»0 
.3730nft(lOF:*0 
.3?*PPnOOF*0 

,3?snoonoE»o 
,3?FpnoooF«n 
,3?7onnnoF»o 
,3?PonoooE*o 
.3?‘*poftooF*<> 
.33000flOOF*0 
.33100P00F*O 
.33?oonooE*o 
.333nftft00F*0 
,334pnft0f!F»0 
.3??iPn(iooF4  0 
.33Fo0"00E*0 
.3?7onnoPF*o 
.33P0onooF*o 
.339PO000F40 
.34nnnfl00E40 
.34]nnnooF»n 
.34?ppnooF*o 
.343PonPOF»0 
.344nnnooE»o 
.34*ipn«onF*n 
.34f OO000F*0 
.347n(innoF*n 
.348ne(iooF*o 
.34«pnnooF*o 
.3‘'OnpnooF»o 
,7«il(lonoOE40 
.3‘>?nnoooF*o 
.3'<300nO0F*0 
3c;4nOnOQF*0 
.3«i<>POP00F*0 

.3«;#t(innooE*0 

.3c;7nnnooF*n 

.3*'«noftOOF*p 

.3‘*‘'noottOE*n 

,3AnnpnftOF»o 


.iflf«nfi*oF*o 

.ioPon4iF*o 

.lpSV0S3ftF.*n 

.  1  o^HOnnnE^O 

.10Sf.9S43F40. 

.10S5VlF.f,E*01 

.]054Hfl77E*n‘ 

.1063ehH3F*0 

.r0528*»<»3E*n 

.1  0S1«F16F*0 

.10S0fi7FnF*O. 

.1049‘>033E*ni 

.10489444E40' 

.io*«opnoE*o 

.i04707nSE40 

.10P61840E*0 

.1  045248 lF*n 

.1044350«E*fl 

. ) p434584E«fl 

.10425F97E40 

.  J04I8P20E*I) 

.)04ft79?7E*n 

.1  0398«»95F*0 

.l0390ftnOE*0. 

.in3809?6E*01 

.in37I79?F*0- 

.1  0382F?8E*0 

.103534F1E*0. 

.10344321F401 

.)0335?34F*0 

.1032*>231E*0 

.1031733«F*0 

.1  030flSHSE»n 

.1  n3bnnonE*o 

.i029if.n?F*n 

.102ft337(lE*0 

,10?75?7f.E*0 

.102ft7?9lE«0 

.in2593«4E*0. 

.102514?f.F.*01 

.10243f.BBE*01 

.J0235«41E*0} 

!i0227955F»0* 

,|02?000nE*O 

.J021l9*'«E*n 

.in203H47F*0 

.J0195F9F.F»n 

.1018753FE*0 

.)OlTi30?E4Ol 
..101 8328  7E*0i 
.10155379E*01 
.30147».n7E*0' 
.ioi4onooE*o 
.1013?sr»E»0 
.101253?3F*0. 
.101I62P8E*01 
.10in?05E*01 
.1P104?B8E*OJ 
.100974?5E*0i 
.10090593F.ni 
.10083783F..01 
.I007690SE*01 
.10070POnF*01, 
.10083ni9E»01 
.10055O72F.01 
.1  n048Fr4F»01 
.ion*i740F.ni 
.100345e4E.0l 
.1002T422F401 
.1 On2n?8BE.01 
. 100131 36E*01 
.lono6n42F*Oi 
.OQ990nP0F»00 
.099202212*00 
. 09851 P98F.00 
,997828122*00 
•  .o9714754E»nO 
,oo8*7S10E*0n 
.995808852*00 
.9951*8052*00 
.99**9317E*0n 
.0938*3872*00 
.993200002*00 
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TABLE  2.  (CONTINUED) 


PT.MO,* 
PT.MO,* 
PT.Mft.» 
PT.nn.s 
PT.NO.a 
PT.NO.m 
PT  ,mO,« 
PT, 

PT.nO.s 
PT.K'O,* 
PT.HO,* 
PT.»»0.= 
PT.NO.* 
PT.Nft.a 
PT..Wft,« 
PT.Mfl.x 
PT.K'O.* 
PT.nn.» 
PT.Nn.s 
PT.NO.s 
PT.NO,* 
PT.MO.« 
PT.Nrt.a 
PT.pn,* 
PT.wn.* 
OT.fP.a 
PT.MP.s 
PT.MP.- 
PT.NO.s 
PT.^n.= 
PT.»'rt.= 
PT.PO.s 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO.* 
PT.NO.* 
PT.NO.* 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO,* 
PT.NO.* 
P>T.Nft,a 
PT.NO.* 
PT.NO,* 
PT.NO.* 
PT.NO,* 
PT.NO.* 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO,* 
PT^NO.* 
PT.NO,* 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO,* 
PT.NO.* 
PT.NO,* 
PT.NO.* 
PT.NO.* 
PT.NO.* 
PT ,NO,= 
PT.NO,* 
PT.NO.* 
•  PT.NO.* 
PT.NO.* 
PT.NO,* 
PT .NO,* 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO.* 
PT.NO.* 
PT.NO,* 
PT.NO,* 
PT.NO.* 
PT.NO,* 
PT .NO,* 
PT.NO,* 
PT .NO.* 
PT.NO.* 


03? 

333 

334 
■*3<i 
33«« 
337 
33« 
33N 

340 

341 
34? 

343 

344 
'34b 

346 

347 

348 
344 

350 

3bl 

Ob? 

OSS 

3S4 

OSS 

356 

357 
OSfl 

359 

360 

361 
36? 
■»63 
36* 

365 

366 

367 

368 

369 

370 

371 
37? 

373 

374 

375 

376 

377 

378 

379 

380 

381 
38? 

383 

384 
386 

386 

387 

388 

389 

390 

391 
39? 

393 

394 
396 

396 

397 

398 

399 

400 
*01 
*0? 
*03 
404 
406 

:r7 

408 

409 

410 

411 
417 
413 

trs 


X*  .361fl0n008»01 
X*  ,36?notinoe*oi 

X*  .osooooooF^ri 

x*  .36*onpooe»oi 
X*  .?6soonooF,*oi 
X*  .?66nooooP*oi 
X*  .36700600P»ni 
X*  ,36800000E*01 
X*  .36«00000p*0-l 
X*  .370000008*01 
X*  ,37iooopoe*oi 
X*  .377000008*01 
X*  ,37300000F*01 
X*  .37400000E*01 
X*  .37SooflnoE*01 

X*  .37600000F*01 
Xb  ,377000008*01 
X*  .37H00000E*01 
X*  .37900000 '*01 
X*  .380000007*01 
X*  .381 A0A00E*01 
X*  .38?POOOOE*ni 
X*  ,38300noOE*01 
X*  .384000flOE*01 

X*  .38snooooE*ni 

X*  .XOOOOOOOE*^ 
X*  .38700000E*01 
X*  .388000008*01 
X*  .389000008*01 
X*  ,3oonooooE*oi 
Xa  .39l000fl0E*01 
X*  .39?OOOOOE*01 
xs  .393006008*01 
X*  ,39400n00E*0 
Xa  ,39SOOOOOP*0 
X*  ,39600000E*0 
Xa  ,39700000E*0 
Xa  ,3O800000E*O 
x*  .399nooooE*n 
x*  ,40onooooE*n 
X*  .401000008*0 
X*  .40?oonooe*o 
X*  .40300<)OOE*01 
x*  .4040Aft00P*0 
X*  .406000908*0 
X*  .406000008*0 
X*  ,40700000E*0 
X*  .40800000F*01 
X*  .4Q900000E*01 
X*  .OIOOOOOOE*!!! 
Xa  .411006008*01 
X*  .41?OOOOOE*01 
Xa  .413ono00E*01 
X*  .414006008*01 
X*  .41SOOOOOE*01 
X*  .41600000F*01 
X*  .417000008*01 
X*  .418006008*01 
Xa  .419000008*01 
X*  .470000008*01 
X*  .471006008*01 
X*  .477000008*01 
X*  .473006008*01 
X*  .474606008*01 
X*  .476000008*61 
X«  .476006008*01 
X*  .477006008*01 
Xa  .4?8n0600E*01 
X*  .*79006008*01 
v=  .*3onnnnnP*nt 
»*  ,*31066008*01 
X4  .*37000008*01 
X*  .433006008.01 
Xs  .434000008.01 
X*  .436000008*01 
X*  .436000008*0] 
X*  .600000008*01 
X*  .600000008*01 
X*  ,700000008*61 
X*  .060666006*01 
X*  .900006008*61 
Xa  •.106066068*07 
X*  ,110060008*0? 
X*  .170060008*0? 


F(X)S 

8  (X)  a 

FIX)* 

8(XJ  = 

FIX)* 

8  IX)  a 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

F  |X)* 

FIX)* 

F|X)s 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

8  I  X  )  « 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

8|X)a 

FIX)* 

FIX)* 

FIX)* 

F|X)a 

FIX)* 

FIX)* 

FIX)* 

F  |X)  = 

f  |X>» 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX)* 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

FIX) 

F|X) 

F  JX) 
FIX) 
FIX) 
FIX) 


.99?56)40E*00 

.99197779E*00 

.991?98N4F*00 

.990674748*00 

.990063668*00 

.989436788*00 

.988873308*00 

.9S8?1789iF.00 

.987605738*00 

.987000008*00 

.986396088*00 

.V86796?8I?.|)0 

.986196178*00 

.984607718*00 

.984010778*00 

.9834?101E*00 

.987836148*00 

.OH7?b?87F*00 

.981674478*00 

.981100008*00 

.OHnS?970F.*00 

.079963108*00 

.979390688*00 

.97B«3F9lf*00 

.978?ftO?SE*00 

.977773178*00 

.977167138*00 

.976611678*00 

.976066008*60 

.976600008*00 

.976943048*00 

. 974386778*00 

.973a?6768*00 

.973?7?868*00 

.97?71H74F*00 

.O7?16713E*h0 

.971610738*00 

.97l074?6e*00 

.970534458*00 

.97O0n000F*00 

.969471438*00 

.96094841)8*00 

.968430198*00 

,967916878*00 

.967407?9F*00 

.«T6690 1138*00 

.966397868*00 

.965096938*00 

.96S397P?E*0n 

.964900008*00 

.r96440306E*00  ■ 

.963907078*00 

,96341??8F.*00 

.96?91888F*00 

.9674?710E*00 

.961937168*00 

.961449?6E*00 

.960963618*00 

.960400468*00 

.960000008*00 

.9696??368*00 

.959047318*00 

.960574518*00 

.960103678*00 

.957634378*00 

.957166?7F*00 

.956699158*00 

.956237678*00 

,956766368*00 

.96530POOF*00_ 

.964033398*00 

.954366098*00 

.953901008*00 

.95343O?4E*00 

.9529731?C*00 

.067517)58*00 

.9?690n008*00 

.096300008*00 

.072700008*00 

.053000008*00 

.037900008*00 

.O?4?000n8*00 

.01196970^*00 

.60095600F*00 
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TABLE  2.  (CONTINUED) 


PT.NO.s 
PT.MO.« 
PT.>'n.a 
PT.^Jft.m 
PT.MO.a 
PT.NO.a 
PT.MO.a 
PT.MO.a 
Pt.MO.a 
PT.MO.a 
PT.MO.a 
.  PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 
PT.MO.a 

PT.Nn.* 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.i 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 

PT.MO.a 


4lf, 

417 

4l« 

41M 

4?0 

421 

42? 

423 

42* 

4?«1 

4?6 

4?7 

42a 

4?9 

430 

431 
43? 

433 

434 

435 

436 

437 
439 

439 

440 

441 
44? 
4*3 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 
469 

460 

461 
46? 

463 

464 

465 

466 

467 

468 

469 

470 

471 
47? 
*73 

474 

475 
*76 
477 
*78 
*79 

490 
*81 
*82 
483 
48* 
*85 

486 

487 

488 

489 
*90 

491 
♦  <>2 
493 
49* 
4V5 
496 
*97 


)ta 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
•  X  = 
Xa 
Xa 
Xa 
Xa 
Xa 
Xa 
xa 
Xa 
Xa 
Xa 
X  = 
Xa 
xa 
X  = 
Xa 
Xa 
Xa 
Xa 
Xa 

x  = 
xa 
Xa 
Xa 
xa 

Xa 

Xa 

Xa 

Xa 

'Xa 

Xa 

X? 

Xa 

Xa 

xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

Xa 

xa 

Xa 

Xa 

Xa 


.13000000^4.6?  F(X)B 

•  .l*flfl0000(-40?  Fix>  = 
. 15o0fl00OP40?  Fix)a 
.160ft0pfl0F40?  FlXia 
.17P0OO00F4O?  FfXia 
.18nO0000F*0?  FlX)a 

.1O000000F40?  F(xia 
.?nnooonnF*o?  Fjxia 

.?! nnO000F*O?  FlXia 
.??a0nO00f 40?  FfX)a 
.?30POOOOF40?  F(X|a 
.?*OOOOOOF40?  FlXia 
.?500(1000F4n?  FiXja 
.?6000000F40?  F(X)a 
.?7noonooE4n?  Fixia 
.?«nonflooe40?  F(xia 
.?90nn000e40?  Fjxia 
.30000000^40?  FIX)B 
.31000000^40?  FlXia 
.3?OOOnoOF40?  FlXl* 
.3.3noonOOF40?  Fixia 
.34000000F40?  FlXl* 
.3S0O0O00e40?  Fixia 
.36n0fl000E40?  FiX)a 
.37000000E40?  FlXia 
.38onoonoE4a?  Fixia 
.39nnOOOOE40?  Fixia 
.40000000F4O?  Fixia 
.*1  onooonf:40?  Fixia 

.*?000O0nE4O?  Fixia 
.43onononF40?  Fixi® 
.44000000F:40?  F|X|a 

.*5nonoooF*n?  fixi= 
.46000000E*Q?  Fixia 
.*7000fl0flF*0?  FIX)* 
.*flOOOOOOE40?  FlXia 
.*9000000E*0?  FlXia 
.5flOOOnOOE*fl?  FlXia 
,5l(loonOOF40?  FIX)  a 
.5?OOOflOOE*0?  FlX)a 
.53n0OO00E*O?  Ftxi» 
.5*oonnooE*n2  Fixia 
.55O0flO0OE*0?  FlX)a 
,560nflOOOE*0?  Fixia 
.57n0Ofl00F*O2  FlX)a 
.5«00fl000f*0?  FlXia 
,59000000F.40?  FlX)a 
.60000000E*0?  FlX)a 
.6ioonoooE*o?  FlXl» 
.6?oonnOQE*n?  Fixia 
.630Ofl00nE*O?  Fixia 
.6*0O0n00F40?  Fixia 
.65o(1OO00E40?  F|X)a 
«66nonnooE«02  Fixi® 
,670Pfln00F40?  FlX)» 
,6POOOOOOE40?  Fix)a 
.69OnOOO0E40?  Fixia 
.7On00000E4O?  FlX)a 
,7ionnnflOE*o?  Fixia 
.7?0PnnOOE4O?  Fixia 
.7300fiO0OE4n?  F|X)a 
.7*00fln00E*0?  FlX)a 
,75nfln00OE40?  F|X)a 
.760onoOOE40?  FlX)a 
.77000000E*0?  FlXia 
.78flO0000E*O?  FlXia 
.79nflOO00E*0?  F1X)3 

.«OnoOOOOF40?  FlXl* 
.81000000E40?  FlX)a 
.fl?OOOnOOE40?  FlX)a 

.P3ooononE*o?  Fixia 
,A*noonooe*o?  Fixia 
.«5oonnooE4n?  Fixi* 

.86000000F40?  Fixia 
,«700n000E*fl?  F1XI» 
.88onnnooF*o?  F|xia 
.fl<lflOOnOOE*0?  Fixia 
.90000nOOE*0?  Fixia 
.9) OOnOOOE402  Fixia 

.9?oonnooF40?  Fixia 
.«3oonnooE*o?  Fixia 

,94000000e40?  Fixia 


.79104803F*00 
.7A213?9*F*00 
.77*09fl34F4n0 
.766H31P7F*0D 
.760?21 lflF*nO 
.75415388E*00 
.74851761E*00 
.74320onOF400 
.73H10636F*00 
.733?1?67E*00 
.7?«5i ?59F400 
.7?399977F*00 
.719667?7F*00 
.7l551O53E*00 
.7115214?E400 
.70769*?OE400 
.7040?250E*00 
.70050O00E4n0 
.6971?006F*00 
.693874«1E*00 
.69075650F400 
.6ft775678F4O0 
.684H6769E400 
.68?08n9F*00 
.6793k9?3e*00 
.6767K374E*0e 
.6742566 8F400 
.671«no00E400 
.6694n67)F*00 
.66707410E400 
.66480n5?F400 
.66?584??F400 
.6604?386F400 
.65H31749F*00 
.65626357F400 
.65426044E*00  ■ 
.6S230647F*00 
.650400POE*00 
.64B53M30F*00 
.64672?30E*00 
.64494684E*00 
.643?1075E*00 
.6415n8TF400 
.63984B04E*00 
.638217)  (1F*00 
.636616P9E*00 
.635045?4E*00 
.63350nonF*00 
.63197938E*00 
.630*8309E400 
.6?90n??F400 
.62756388E400 
.6?61*n6E*O0 
.624743)4E*00 
.623369V?E*00 
.62202)602*00 
.620698?6F*00 
.61940onOF*00 
46)812670F*00 
.6)667735X400 
.6l565n76F400 
.61444571E*00 
.61326n99E*00 
.61?09639E*00 
.61094771E400 
.6098167?F*0(1 
.60670)22E*00 
.607600POE*00 
.60651204E*00 
.60543710E*00 
.604375)3EtOO 
.603326n7E*00 
.6022e987F*00 
.60) 2664aF*00  . 
.6n02b585E*00 
.599?5793E4n0 
.59fl27?66F*00 
.5973000nE*00 
.5963398 3E*00 
.59539)83F*C0 
. 5944556 ?e*00 
.5935308)E*00 
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TABLE  2.  (CONCLUDED) 


I 


PT.NO.s 

4PB 

Ifs 

.pspooftooe^o? 

FIX) 

P7.9n*= 

499 

yx 

.9#.0O0ftn0F»0? 

FIX) 

PT.Wrt.s 

600 

yc 

.O7000000F*0? 

FIX) 

PT.NO.a 

601 

y« 

.OH000ft00E*0? 

FIX) 

P7.N0.'» 

60? 

Xx 

.99onooooF*o? 

.  FIX) 

PT.NO." 

603 

x« 

.10flOn000E*O3 

FIX) 

..S9?617niF*00 

.«191713(»BE*00 

.•;90S?098F»00 

.*;«99379ftE*00 

.SB9Q89B3E*00 

.!>S8?0000^*00 


IXXX08L 


been  calculated  by  Svehla  |2|  for  O;  and  NO.  However,  a  literature  search  did  not  reveal 
viscosities  of  Oj  and  NO:  at  elevated  temperatures. 


A  method  for  the  calculation  of  the  viscosity  of  gases  is  given  by  Bird,  et  al.  [3] 


e/k 

=  0.77T 

c 

(11) 

/T 

a  = 

2-44  (^) 

(12) 

\  / 

0  = 

~  Va 

0.841  (V^) 

c 

(13) 

where  these  functions  are  utilized  in  the  Lennard-Jones  potential 


(J>(r)  ®  46 


(14) 


For  O3 

=  268°K  Pc  =  67  ATM 

V  =89.4  cmVgro  mole  MW  =  48.000 
c 

For  NO2 

=  431.0  °K  =  100  ATM 

G  C 

MW  =  46.008 


Utilizing  (II)  and  (13)  for  O3  and  (II)  and  (12)  for  NO:  produced  the  following 
potential  functions  for  O3  and  NO:; 


i 
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GAS 

e/k 

0 

MW 

NO2 

331.8 

3.97 

46.008 

O3 

206.4 

3.76 

48.000 

The  viscosity  for  each  of  these  gases  is  then  calculated  using: 
y  =  2.6693  x  lO"^  {iM)__T}  h 

Utilizing  the  potential  functions  given  in  Table  2,  Equation  (15),  and  the  program 
developed  in  the  previous  section  to  calculate  Cl^  gives  the  viscosities  of  O3  and  NO2. 


A  program  was  written  to  accomplish  this  called  MUCALC  and  a  listing  of  this 
program  is  given  in  Appendix  B. 

Table  3  gives  the  results  of  these  calculations  for  NO2  as  a  function  of  temperature 
while  similar  results  for  O3  are  given  in  Table  4.  Finally  results  for  the  other  gases  of  interest 
that  have  been  taken  from  the  literature  are  presented  in  Table  5.  These  gases  are  N2,  NO,  and 
O2.  These  results  are  given  over  the  narrower  range  of  temperature  of  interest  for  these 
reactions. 


IV.  COMPARISON  OF  RESULTS  OF  THIS  METHOD  WITH  OTHER 
METHODS  FOR  CALCULATING  MOLECULAR  VISCOSITIES. 

In  order  to  determine  the  accuracy  of  the  technique  utilized  to  calculate  molecular 
viscosities  that  were  presented  in  the  last  section,  a  calculation  was  made  for  CO2,  results  for 
which  can  be  found  in  the  literature.  Hence,  the  properties  for  CO2  were  introduced  into  the 
program  presented  in  the  preceding  section  and  the  molecular  viscosity  for  CO2  was 
calculated.  The  program  listing  for  this  is  presented  in  Appendix  C.  Note  that  this  program  is 
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TABLE  3.  MOLECULAR  VISCOSITY  OF  NO2  AS  A  FUNCTION  OF 
TEMPERATURE 


Ts 

?00 .0 

OFf> 

K 

y  V|i= 

T  = 

3on,o 

OFF 

k 

y  vii= 

T= 

400.0 

nc-R 

tc 

y».Mj= 

T  = 

soo.o 

OFF 

y '•11  = 

T= 

bOO.O 

OFF 

K 

yM»)= 

T= 

700.0 

OFF 

K 

y-iia 

.  T= _ 

Hoa.n 

0F6 

K 

XMlls 

T= 

900.0 

OEF 

K 

HNIIa 

T  = 

1000. 0 

off 

K 

y  Mila 

T  = 

1100.0 

OFF 

y 

XM|I  = 

T= . . 

1200.0- 

OFF 

K 

XPlia 

T= 

1300.0 

npp 

K 

XMIJS 

T= 

1400.0 
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TABLE  4.  MOLECULAR  VISCOSITY  OF  O3  AS  A  FUNCTION  OF 
TEMPERATURE 
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TABLE  5.  MOLECULAR  VISCOSITIES  FOR  GASES  INVOLVED 
IN  THE  N0+03-^N02  +  O2  REACTION  IN  A  N2 
CARRIER  GAS. 


identical  to  the  program  presented  in  Appendix  B  except  for  the  input  conditions  which  are 
shown  highlighted. 

The  comparison  of  the  results  from  this  method  and  that  of  Svehla  are  shown  in  Table  6 
for  CO2.  The  differences  between  the  results  range  from  1.4  percent  at  200  degrees  K  and 
increase  monatonically  to  3.3  percent  at  5000  degrees  K.  Hence  the  conclusion  is  that  the 
calculational  method  results  in  small  deviations  from  accepted  results  for  gases  for  which  the 
more  vigorous  treatment  has  been  exercised. 


V.  INTERPOLATION  FOR  m.  =  Mi(T) 

Having  now  established  an  interpolation  routine  for  specie  viscosities  and  having 
calculated  specie  viscosities  for  NO2  and  O3  plus  estimating  the  general  magnitude  of  the  error 
involved  it  was  necessary  to  establish  that  the  interpolation  procedure  was  working  correctly 
for  individual  species.  This  was  accomplished  by  programming  MUSPEC  which  is  given  in 
Appendix  D.  This  program  utilizes  the  viscosities  of  N2,  O2,  NO,  NO2,  and  O3  as  input  and 
calculates  the  viscosities  of  these  species  at  various  temperatures  between  500  degrees  K  and 
600  degrees  K. 

The  results  of  these  calculations  are  shown  in  Table  7.  These  results  can  be  compared 
with  the  data  presented  in  Table  5.  This  comparison  shows  that  the  calculations  are 
consistent  and  reasonable.  Hence,  the  procedure  for  the  calculation  at  the  individual 
viscosities  has  been  verified  before  proceeding  to  the  mixture  of  these  gases  given  in  the  next 
section. 


VI.  CALCULATION  OF  THE  MOLECULAR  VISCOSITY  OF  A 
MIXTURE  OF  GASES 


Having  established  the  calculational  procedure  for  individual  viscosities  of  gas  species, 
the  remaining  task  is  to  calculate  same  for  the  mixture  resulting  from  the  chemical  reaction. 
The  viscosity  of  the  mixture  of  gases  is  not  simply  a  mole  fraction  average  of  the  individual 
viscosities  but  depends  on  the  individual  species  in  a  more  complex  manner.  The 
mathematical  model  due  to  Wilke  is  given  [3]  as  follows 
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=  I 
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(16) 


TABLE  6.  COMPARISON  OF  THE  VISCOSITY  CALCULATED  BY  THE  PRESENT 
METHOD  AND  THAT  OF  SVEHLA. 


MOLECULAR  VISCOSITY  OF  CO2  AS  A  FUNCTION  OF  TEMPERATURE 
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where 


Equation  (16)  has  been  shown  to  reproduce  measured  values  of  iimx  to  within  2  percent 
for  certain  gases.  The  dependence  of  Hmx  on  composition  is  extremely  non-linear  for  certain 
gas  mixtures,  however. 

Nevertheless,  it  is  the  best  available  technique  and  will  be  utilized  here. 

A  test  program  was  written  to  compute  the  viscosities  of  various  mixtures  of  gases.  The 
listing  of  this  program  is  given  in  Appendix  E.  This  program  utilizes  subroutine  LAM  VISC 
which  performs  the  mixture  calculations  for  the  laminar  viscosities  according  to  the  technique 
of  Wilke  given  above.  It  also  utilizes  the  subroutine  M  U  SPEC  which  is  a  variation  of  Program 
M  US  PEC  discussed  in  an  earlier  section. 

The  results  utilizing  these  methods  are  given  in  for  various  mixtures  of  N:,  O2, 

NO,  NO2,  and  O3.  Comparing  table  8  and  Table  5  will  convince  one  that  the  results  are 
reasonable. 

VII.  CONCLUSIONS 

A  computer  code  has  been  generated  to  determine  laminar  viscosities  of  gases  as  a 
function  of  temperature.  This  method  is  approximate  but  it  provides  results  useful  in  making 
engineering  analyses.  The  specific  gases  for  which  laminar  viscosities  were  determined 
included  NO2  and  O3,  but  the  method  is  applicable  for  other  gases. 
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TA^LE  8.  Mmix  =  Mmix  (T)  FOR  VARIOUS  MIXTURES  OF  N2, 02.  NO,  NO2.  and  O) 
(MOLE  FRACTIONS) 
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APPENDIX  A 

PROGRAM  LISTING  FOR  TEST  -  A  ROUTINE  TO  CHECK  THE 
ACCURACY  OF  THE  CUBIC  SPLINE  INTERPOLATION  ROUTINE 
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APPENDIX  B 

PROGRAM  LISTING  FOR  MUCALC  -  A  ROUTINE  TO 
CALCULATE  THE  MOLECULAR  VISCOSITY  OF  NO2 
AND  O3  FROM  200  DEGREES  K  -  5000  DEGREES  K 
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APPENDIX  C 

PROGRAM  LISTING  FOR  MUCALC  SPECIALIZED  FOR 
CALCULATION  OF  THE  MOLECULAR  VISCOSITY  OF  CO2 
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APPENDIX  D 

PROGRAM  LISITNG  FOR  MUSPEC  -  A  ROUTINE  WHICH 
CALCULATES  n  =  Mi(T)  FOR  N2,  O2,  NO2,  NO,  AND 
O3  AT  TEMPERATURES  BETWEEN  200  DEGREES  K  AND  1000 

DEGREES  K 
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SUBROUTINE  CALCCF  -  IDENTICAL  TO  CALCCF  IN  APPENDIX  A 


APPENDIX  E 

PROGRAM  LISTING  FOR  LVSCTST  -  A  ROUtiNE  TO 
CALCULATE  mmix  =  Mmix(T)  FOR  MIXTURES 
OF  Nj,  O2,  NO,  NO2,  AND  O3 
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SYMBOLS 

Viscosity  potential  function  of  (IfT/t) 

Characteristic  energy  of  interaction  between  molecules,  erg/ molecule 
Intermolecular  distance,  cm 

O 

Collision  diameter  of  a  molecule,  A 

Critical  temperature,  degrees  K 

Critical  pressure,  atm 

Critical  volume,  gm/gm-mole 

Boltzman  constant,  I.380S  erg/ molecule-degrees  K 

Lennard-Jones  potential.  Equation  (14) 

Molecular  weight,  i***  specie,  gm/gm-mole 
Molecular  viscosity,  i"*  specie,  poise 
Mole  fraction,  i"'  specie,  dimensionless 
Viscosity  weighting  function.  Equation  (2) 
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